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Abstract 

By considering the x-dependence of Tr"*", vr", , K~ , A, A, p, p hadron productions in charged 
lepton semi-inclusive deep inelastic scattering off nuclear target (using Fe as an example) and 
deuteron D target, we find that (A'^/A'^)/(A'^/A^) and {p^/p^)/{p'^/p^) are ideal to figure out 
the nuclear sea content, which is predicted to be different by different models accounting for the 
nuclear EMC effect. 
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I. INTRODUCTION 



In 1982, the European Muon Collaboration (EMC) at CERN found that the structure 
function ratio of bound nucleon to free nucleon, in the form of F.^ {x , Q"^) / F^f {x , Q^) , is 
not consistent with the expectation by assuming that a nuclei is composed by almost free 
nucleons with Fermi motion correction taken into account [l, 2], and such phenomenon was 
confirmed by E139 collaboration at SLAC [3]. This discovery, which is called the nuclear 
EMC effect, has received extensive attention by the nuclear and hadronic physics society. 



and the rescaling model 



fly, 
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Many nuclear models, such as the pion excess model [4, [s], the quark cluster model 
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12| et ai, have been proposed to explain the data, and 



all these models can qualitatively describe the data in the mediate x region. The inclusive 
deep inelastic scattering (DIS) data are expressed as F^{x,Q'^)/F^{x,Q^), which can be 
written in the naive parton model as: 



F^^{x,Q^) _ S,ef fe(x,g^A) + g,(a;,Q^A)] 
F.^(x,Q2) 



(1) 



where Cj denotes the charge of the partons with flavor i, and q{x, Q"^) is the parton distri- 
bution function of a nucleon. Fig. [1] shows the F^{x,Q^)/F.p{x,Q'^) results of the cluster 
model, the pion excess model and the rescaling model respectively at = 5 GeV^ in the 
mediate x region. All these models, as can be seen from Fig. [1], predict similar behavior of 
Q^) I F.p{x, Q^) at mediate x region. However, the sea quark of the nuclei is differently 
described by the three models. In the cluster model, all sea quarks are enhanced. In the 
pion excess model, the sea quarks u and d are enhanced while the other quarks are reduced. 
However, in the rescaling model, all sea quarks are reduced in the nuclei compared with 
those in the free nucleon (Fig. [2]). 

The Fermilab experiment 772 13| measured the dimuon yield in Drell-Yan process in- 
duced by 800 GeV proton off various nucleus and compared the data with the theoretical 
predictions of the three models accounting for the EMC effect [l^. The data were ex- 
plained in Ref. 13| to favor a conclusion that the sea quark in the nuclei is not enhanced, 
by neglecting the energy loss effect of the incident quark, which is not precisely determined 
yet 1% llSl, llGj. To avoid the uncertainties concerning the sea quarks in the nuclei by the 
dimuon yield in Drell-Yan process solely, the sea content can be also measured in other ex- 
periments. The purpose of this work is to show that the semi-inclusive hadron productions 
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FIG. 1: Results of F^{x)/F^{x) in three models at = 5 Qev2. The sohd, dashed and dot- 
ted curves are the results of the cluster model, the rescaling model and the pion excess model 
respectively. The target nuclei assumed here is Fe. 



in charged lepton deep inelastic scattering are sensitive to the sea quark content of the nu- 
clei. The pion, kaon, proton and antiproton productions in the charged lepton semi-inclusive 
deep inelastic scattering off nuclei have been checked by the HERMES collaboration jl^. 
However, while the HERMES data, which are expressed as multiplicity ratios and also with 
2;-dependence, are convenient to study the modification of the fragmentation functions in 
nuclear environment they are not ideal to provide much information about the sea 
quarks of the nuclei. We will show that the A/A production, especially the x dependent 
behavior, is ideal to distinguish between different predictions on the sea content of the nuclei 
in the different models of the EMC effect. 



II. NUCLEAR MODELS AND THE SEA QUARK DISTRIBUTIONS 



The semi-inclusive hadron productions in charged lepton deep inelastic scattering can be 
related with the quark distribution functions as: 



dxdydz 



;i + (1 - yf) J2 e- [q^{x, Q')DI{z, Q') + Q')DI{z, Q^)] , (2) 
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FIG. 2: u and s sea quark behaviors of various models. The solid, dashed and dotted curves 
correspond to the predictions of the cluster model, the rescaling model and the pion excess model 
respectively. The parton distribution is calculated at = 5 GeV^. The target nuclei assumed 
here is Fe. 

where Q^) is the parton distribution for quarks with flavor i, and D^.i^ZjQ"^) is the 
fragmentation function of quark to hadron h. The formula is also applicable to the nuclei, 
with the parton distributions and fragmentation functions replaced by Q^)— Q^, A) 
and D{z, Q'^)—>-D{z, Q^, A) respectively. The inclusive production itself can not offer enough 
information about the sea quark enhancement, while the ratio 

da^Jdx J' dzZ.eKqfix, Q')D^S^, Q\ A) + gf (x, Q^)Dl{z, Q\ A)) 



da^n/dx dz^Z^e^qFix, Q')D^S^, Q^) + gf (x, Q^)Dl {z, Q^)) 



(3) 



is useful to reveal the difference between the sea quark behavior in the nuclei and that in 
the nucleon. 

In the following, we will consider the ratios {da\/ dx) / {da^/ dx) for various hadrons. We 
use the pion excess model, the rescaling model and the pion excess model to calculate the 
sea quark content of the nuclei. For the sea quark distributions in the deuteron D, we use 



the result offered by the model itself 



a 



for the cluster model, and for the other two models 



we adopt the CTEQ5L parametrization [18] of parton distributions for free nucleons by 
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considering the isospin symmetry between proton and neutron. 

In the pion excess model, the quark distribution in the nuclei is modified by the extra 



pions caused by the interaction between the nucleons in nuclei {4]. The quark distribution 
of nuclei is: 

W = ^f.^ym-) + ^/^(i/)gf (-), (4) 
Jx y y Jx y y 

in which qf{x) and qf{x) are the parton distributions in the free pion and in the free nucleon 
respectively, and fn{y) is the probability to find extra pions in the nuclei j^. For simplicity. 



19| . in which the proton is supposed to be 



we adopt the parametrization in a toy model 

partially in the nucleon-pion subsystem state and the parton distributions in the nucleon 

and in the pion are assumed to be the same as those in the free nucleon and in the free pion. 

□ 

Thus, the excess pion and the nucleon probabilities per nucleon are given as [19|: 

fUz) = {1 - {n,mz - 1) + f,^{l - z), (6) 



where 71^ = 0.22, a = 1 and 6 = 3. The CTEQ5L [181] parametrization of the parton 







distribution of the nucleon and MRS [20| parametrization of the parton distribution of the 
pion are adopted to obtain F^{x) and the parton distribution of the nuclei. 

In the quark cluster model, six or more quark cluster is supposed to exist in the nuclei to 
account for the EMC effect. For the sake of simplicity, only six quark cluster is considered 
here. q{x) in a six-quark cluster can not be measured directly from experiment, but Carlson 
and Havens [2] estimated quark distribution q(x) per nucleon in the six quark cluster based 
on QCD counting rules: 

v,ix) = N,z^-'/'\l-zY^, (7) 
ue{x) = (iV,,,/4)^(-l)(l-^)l^ (8) 

where = 1.3875 and Nsea = 0.2521 are the coefficients to warrant momentum con- 
servation and the quark number. In the six quark cluster, x, the variable defined by 
(5^/(2Mjvz/), equals to 2z because z is defined by Q'^I{2Mqu) [7]. Therefore, and 
F^{x,Q'^)/Fl^{x,Q'^) can be given as: 

q^ix) = (l-f)q^(x) + fq%x), (9) 



f) + f JnMZ..^ ^ (10) 



where / is the probabihty to find the six-quark cluster in the nuclei and is adjusted to fit 
the inclusive deep inelastic e{fi)-A scattering data. Its value is given as 0.30 by Carlson and 
Havens 

For the rescaling model, the quark in the nuclear medium is considered to have different 
confinement size compared with that of the quark in the free nucleon. Q^) is related 

with Q^) (the parton distribution in the free nucleon) by the relation: 

g^(x,Q2) = g^(x,e(gW, (11) 



where ^(Q^) varies with A and Q^. ^ equals to 1.83 while = 5 GeV^ and A=56 (Fe). 
For q^{x, Q"^), the parton distribution per nucleon, CTEQ5L parametrization [18!] is adopted. 

Given the above analysis, F^(x,Q^)/ F^(x,Q^) and sea quark enhancement are checked 
in Fig. □ and Fig. H 

III. FRAGMENTATION FUNCTION AND PARTON ENERGY LOSS MODEL 

Due to the non-perturbative nature of the fragmentation process, the fragmentation func- 
tion can not be calculated from first principle, thus models are used to obtain the fragmen- 
tation function. Experimentally, the process of e"*" + e~ h + X can offer much information 



about the fragmentation 21 1. 



Based on the experimental data and theoretical analysis, Kretzer 22| gave a parametriza- 
tion of the fragmentation in the form: 

D^S^,Ql) = Nx^{l~xf, (12) 

where A^, a and (3 are the constants chosen to fit the experimental data. For specified hadron, 
a, which determines the low z region behavior of the fragmentation function, is the same for 
all light fiavor quarks while /3, related with the high z region behavior of the fragmentation 
function, is different for various quarks. Therefore, in low z region, all quark fragmentation 
functions have the same shape, and in the large z region, the favored quark fragmentation 
function is larger than the unfavored quark fragmentation function. In addition, because of 
the strange dominance at large x for s in K~ , the fragmentation function such as s — > K~ 
is larger than that of m K~ . 



As for A, A, p, there exists a phenomenological parametrization 23|] of their fragmenta- 



tion functions based on the assumption that the fragmentation function of quark q to hadron 



h is proportional to the q quark distribution in the hadron h: 

D'^iz) (X q\z). (13) 

In general, the fragmentation functions can be written as: 

D^{z) = Cy{z)z^qUz). (14) 
D^siz) = Cs{z)z^q^s{z\ (15) 

where D^gi^z) means unfavored fragmentation function. There are three options for the 
favored and unfavored quark fragmentations: (1). Cy = 1 and Cs = for a = 0; (2). 
Cv = 1 and C5 = 1 for a = 0.5; (3). Cy = 1 and C5 = 3 for a = 1. The parton 
distributions of A and A are essential to get the fragmentation functions of quark to A 
and A. Based on the fact that there is no direct parton distribution of A and A, SU(3) 



symmetry between the proton and the A is adopted to get the parton distribution in A [23|. 
Although there are three models, SU(3) symmetry model, quark-diquark model and pQCD 
based analysis, the difference of the fragmentation functions will not affect the qualitative 
predictions by the fact that only the flavor structure of the parton distributions of the proton 
are different in the three models in the high x region [2J]. 

The common feature of the above parameterizations of fragmentation functions is that 
the favored quark fragmentation is larger than the unfavored quark fragmentation in the 
large z region. Thus the favored quark fragmentation process is able to obtain sea quark 
information of the nuclei if the produced hadron is from the favored fragmentation of sea 
quarks in the nuclei. Then the hadron events of produced particles in the large z region 
can be chosen to get the x dependence of the production ratio, in which the unfavored 
fragmentation contribution from the valence quark can be largely suppressed. 

Fragmentation function in the nuclei is important for producing the hadron production 
from the nuclei. HERMES collaboration has measured the hadron production from the 



nuclei and found that the production is reduced compared with that from the 

n r 

cleon [17[, and many effects such as nuclear absorption i25i|], parton energy loss 
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ree nu- 
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291] have been developed 



gluon bremsstrahlung [28| and partial deconfinement 
to account for the data. In this paper, the parton energy loss model is adopted to get the 
fragmentation function in the nuclei. 

In Refs. {26! . I27I . the modification of the fragmentation function is caused by the inter- 
action between the hard quark and the debris of the nuclei. Given that the original parton 



energy loss model is complicated to apply in the real process, an effective model suggested 



in Ref. SOj is used in Ref. [311] to get the modification of the fragmentation function. In 
effective parton energy loss model, modified fragmentation function is expressed in the form: 

Q\A)= t ""'^ deD{e, u)z*D^^{z\ Q^), (16) 

^0 

where = fi — e. Eh is the measured hadron energy and e is the energy loss of the hard 
quark going through the nuclei, z* is the rescaled momentum fraction caused by the quark 
energy shift in presence of QCD medium: 

Di^tjQ"^), the probability for a quark with energy E = u to lose energy e, is parameterized 



by Arleo [3^: 



where fi, a are two parameters with a = 0.73 and fi = —1.5, as the energy of the quark, 
which has absorbed the virtual photon, is much higher than the energy loss when it passes 
through the nuclei environment. And uJc is the relevant scale of the typical gluon energy 
and denotes the energy loss scale of the hard quark, 

UJc = \qL\ (19) 

Although q is not precisely determined yet, in the next section we will show that the result 
is not sensitive on q when is large enough. Here we set g = 0.72 GeV/fm^ and L = 3/4i?, 



where R is the nuclear radius 



3l|. 



From another point of view, the modification of the fragmentation function in the parton 
energy loss model is from assumption. Its conformation with the theoretical framework of 
factorization and renormalization is not fully justified, as the definition of the fragmentation 
functions are vacuum matrix elements with no relation to the target material. As we will find, 
the qualitative conclusion of our paper on the ratio (A'^/A'^) / (A-^/A^) will not be influenced 
by including the modification of the fragmentation function in the nuclear environment. In 
order to get rid of the nuclear absorption or the energy loss process, it would be fine to 
consider hadron production at larger energy (say larger than 20 GeV), hence at higher 
virtuality, where we know from the present HERMES data that these two effects prove 
negligible. 
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FIG. 3: (Color online)The production ratios of tt^,tt^ at = 5 GeV^ in the region 0.7 < z < 0.9 
and 0.3 < z < 0.5. The left two figures are the results of 7r+ and the right are vr^. The upper 
two figures are in low z region and the lower two are in high z region. The solid, dashed and 
dotted curves denote the results of the cluster model, the rescaling model and the pion excess 
model respectively. The thick curves denote the production with all quark fragmentation and the 
thin curves correspond to the results with only favored quark fragmentation. The left two are 
TT~^^ /tt'^^ , and the right two are tt^^/tt^^ The target nuclei assumed here is Fe. 

IV. RESULTS 

Fig. [3] and Fig. H] present the results of TT^{ud), 7T^{ud) and respectively. They show 
that there is no large difference between various model predictions in the large z region and 
in the low z region, no matter by including the favored quark fragmentation process only or 
by including all favored and unfavored fragmentation processes. The reason is that ^'^{ud), 
-K^iud) and are contributed mainly by the favored fragmentation processes u 7r+, 
d ^ TT^ and u because that the valence quarks u and d are predominant over the sea 

quarks in the mediate x region and that production of those hadrons are dominated by the 
behavior of the valence quark in the nuclei (Fig. [5]). 

When focused on K~ , the result with only favored quark fragmentation functions is 
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FIG. 4: (Color online)The production ratio of at = 5 GeV^ in the region 0.7 < z < 0.9 
and 0.3 < z < 0.5. The left is the result in the low z region and the right is the result in the 
high z region. The solid, dashed and dotted curves are the predictions of the cluster model, the 
rescaling model and the pion excess model respectively. The thick curves denote the production 
with all quark fragmentation and the thin curves correspond to the results when only favored quark 
fragmentation. The target nuclei assumed here is Fe. 

different from that when all fragmentation processes are considered (Fig. [6]), which indicates 
that the unfavored quark fragmentation function, u and d to K~ , can not be neglected both 
at low and high z region for the predominance of the valence quarks in the x region we 
considered. In the high z region, due to parton energy loss, the K~ is largely suppressed 
and we can hardly see any difference in the three model predictions with all fragmentation 
processes being considered. 

A and A production ratios in different z regions are examined and the results are plotted 
in Fig. [7| and Fig. [HI For the same reason as K'^, three models predict almost the same 
x-dependence of A^/A^. From Fig. [HI difference are generated among various model pre- 
dictions on the x-dependence of A^/A^, and the these difference are not sensitive to the 



three options of the fragmentation functions 



. The reason is that the dominant produc- 



tion of A is through the favored fragmentation of anti-quarks inside the targets, so that the 
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FIG. 5: The ratio u^{x) /u^ {x) of the three models at Q^=b GeV^. The sohd, dashed and dotted 
curves are the predictions of the cluster model, the rescaling model and the pion excess model 
respectively. The target nuclei assumed here is Fe. 

x-dependence of the production ratio is sensitive to the sea quark behaviors of the nuclei. 
But in the large z region, due to the parton energy loss effect, the large difference between 
three models do not manifest themselves significantly as expected. The reason is, at large z 
region, the fragmentation function is largely modified by the parton energy loss. And such 
a phenomenon also happens on A^/A^. Fig. [8] shows that A'^/A^ is not an ideal variable to 
figure out the nuclear sea quark content by including the modification of the fragmentation 
function in the nuclear environment. Fortunately, largely difference between three models 
appears for the quantity (A'^/A'^)/(A'^/A^), which is more accessible in experiment than 
A^/A^. When the integral upper limit is fixed with unity and the lower limit varies from 0.4 
to 0.7, (A^/A^)/(A^/A^) is still model dependent (Fig. [9]). Thus it is plausible to conclude 
that (K^ / K^) / (K^ / K^) can offer information about the sea content of the nuclei. 

Similar to A'^/A^, the antiproton and proton production ratio {p^)/{p^) can not offer 
much information while p^/p^/p^/p^ do generate large difference with different nuclear 
model (Figs. [T0lfT2|) . So, the ratio {p^ / p^) / {p^ / p"^) is another choice to check the sea 
content of the nuclei in experiment. Attention should be paid to extract possible background 
contribution as a large number of protons and antiprotons might be produced from the decays 
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FIG. 6: (Color online)The production ratio of K~ at = 5 GeV^. The solid, dashed and 
dotted curves are the predictions of the cluster model, the rescaling model and the pion excess 
model respectively. The left figure is the result in the region 0.3 < z < 0.5 and the right is the 
result in the region 0.7 < z < 0.9. The thick curves denote the total results including all quark 
fragmentation processes and the thin curves correspond to the results when only favored quark 
fragmentation processes. The target nuclei assumed here is Fe. 

of other baryons. 

g is a sensitive parameter that could largely affect the modification of fragmenta- 
tion function and is not determined clearly yet. We also calculate A^/A^, A^/A^ and 
(A^/A^)/(A^/A^) when q = 0.25 GeV/fm^ (Figs. ^ME). From the fi gures we can con- 
clude that A^/A^ and A^/A^ are largely affected by different q, while (A^/A^)/(A^/A^) 
is almost q independent. And {p^ / p^) / {p^ / p^) also has such properties(Figs. [TT]fT2l) . So, 
the ratio (A^/A^)/(A^/A^) and {p^ / p^) / {p^ / p^) are not sensitive to q. 

Besides, we should mention that the fixed order calculation used here is not appropriate to 
describe hadron production at high z where large logarithms need to be resumed. Therefore 
we should consider the results here as qualitative predictions rather than quantitative ones. 
For more convinced quantitative predictions, we would need better constrained fits of the 
nuclear parton distributions, rather than the earlier EMC model results adopted in this 



12 




FIG. 7: (Color online) The ratio A^/A^ at = 5 GeV^. The sohd, dashed and dotted curves are 
the predictions of the cluster model, the rescaling model and the pion excess model respectively. 
The left figure is the result in the region 0.3 < z < 0.5 and the right is the result in the region 
0.7 < z < 0.9, with the thin, normal and thick curves corresponding to the three options of the 
unfavored fragmentation (1), (2), (3) respectively. The target nuclei assumed here is Fe. 

paper. 



V. SUMMARY 



In this paper, we adopted three models of the nuclear EMC effect: the cluster model, 
the rescaling model and the pion excess model, to calculate their predictions on the hadron 
production ratio in charged lepton semi-inclusive deep inelastic scattering off nuclei in the 
large z region. Our purpose is to find hadrons which are produced mainly from the sea quarks 
of nucleus, so that we can distinguish between different predictions on the sea content of 
the nuclei. For completeness, we considered the production ratios of n^'^/n^^, n^^/n^^, 
X+Ajx+D^ K-^/K~^, A^/A^, A^/A^,(A^/A^)/(A^/A^), {p^/p^)/{p^/p^) and found 
that the ratios of {K"^ / K^) / {hP / IsP) and are ideal to figure out the sea 

content of the nuclei. 
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FIG. 8: (Color online)The production ratio of K^/k^ at = 5 GeV^. The solid, dashed and 
dotted curves are the predictions of the cluster model, the rescaling model and the pion excess 
model respectively. The left figure is the result in the region 0.3 < z < 0.5, the middle is the result 
in the region 0.4 < z < 0.9 and the right is the result in the region 0.7 < z < 0.9. The thin, 
normal and thick curves correspond to the three options of the unfavored fragmentation (1), (2), 
(3) respectively. The target nuclei assumed here is Fe. 

More significantly, (K"^ / K"^) / (K^ / K^) and {p^ j "p^^ j {p^ j -p^) are accessible in experiment 
and the behaviors of {hJ^ j bJ^) j {bP j bJ^^ and {p"^ / p^) j {p"^ j p"^) are different for different 
models. According to Fig. [9] and Fig. [12] , we conclude that the various models about the 
EMC effect with different sea behaviors can be distinguished by the future data of the 
x-dependence of (A'^/A'^)/(A'^/A^) and {p"^ / p"^) / {p"^ / p"^) in semi-inclusive deep inelastic 
scattering process. The difference between the pion excess model and cluster model is not 
good enough to be checked out in experiment, but whether the sea quark is enhanced or not 
is clear to be distinguished according to the result. 
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FIG. 9: (Color online) The production ratio of (A^/A^)/(A^/A^) at = 5 QeV^, calculated in 
0.7 < z < 1.0 and 0.4 < z < 1.0. The solid, dashed and dotted curves are the predictions of the 
cluster model, the rescaling model and the pion excess model respectively. The thin, normal and 
thick curves correspond to the three options of the fragmentation (1), (2), (3) respectively. The 
target nuclei assumed here is Fe. 

of Higher Education (China). 



[1] J.J. Aubert, et al., Phys. Lett. B 105,322 (1982). 

[2] CERN NA2/EMC, J.J. Aubert, et al., Phys. Lett. B 123, 275 (1983). 

[3] E139, R.G. Arnold, et.al., Phys. Rev. Lett 52, 727 (1984). 

[4] C.H. Llewellyn Smith, Phys. Lett. B 128, 107 (1983). 

[5] M. Ericson and A.W. Thomas, Phys. Lett. B 128, 112 (1983). 

[6] R.L Jaffe, Phys. Rev. Lett. 50, 228 (1983). 

[7] C.E Carlson, T.J. Havens, Phys. Rev. Lett. 51, 261 (1983). 

[8] H.J. Pirner and J.P. Vary, Phys. Rev. Lett. 46, 1376-1379 (1981) 

[9] F.E. Close, R.G. Roberts, and G.C. Ross, Phys. Lett. B 129, 346 (1983). 

[10] R.L. Jaffe, et al., Phys. Lett. B 134, 449 (1984). 

15 




FIG. 10: (Color online) The ratio ^^/p^ (right) and p^/p^(left) at = 5 q^^^ rpj^^ ^^j^^^ 
dashed and dotted curves are the predictions of the cluster model, the rescahng model and the 
pion excess model respectively. Thin, normal and thick curves correspond to the three options of 
the unfavored fragmentation (1), (2), (3) respectively. The target nuclei assumed here is Fe. 

[11] F.E. Close, R.L. Jaffe, and R.G. Roberts, and G. G. Ross, Phys. Rev. D 31, 1004 (1985). 

[12] O. Nachtmann and H. J. Pirner, Zeit. Phys. C 21, 277 (1984). 

[13] D.M. Aide, et. al., Phys. Rev. Lett. 64, 2479 (1990). 

[14] M.A. Vasiliev, et. al., Phys. Rev. Lett. 83, 2304 (1999). 

[15] M.B. Johnson, et. al., Phys. Rev. Lett. 86, 4483 (2001). 

[16] G.T. Garvey and J.C. Peng, Phys. Rev. Lett. 90, 092302 (2003). 

[17] HERMES Collaboration, A. Airapetian, et. al., Phys. Lett. B 577, 37 (2003). 

[18] H.L. Lai, et. al., Eur. Phys. J. C 12, 375 (2000). 

[19] E.L. Berger, F. Coester, and R.B Wiringa, Phys. Rev. D 29, 398 (1983). 

[20] P.J. Sutton, A.D. Martin, R.G. Roberts, and W.J. Stirling, Phys. Rev. D 45, 2349 (1992). 

[21] O. Biebel, P. Nason, and B.R. Webbler, |hep-ph/0109282 

[22] S. Kretzer, Phys. Rev. D 62, 054001 (2000). 

[23] B.-Q. Ma, I. Schmidt, and J.- J. Yang, Phys. Lett. B 574, 35 (2003); 

B.-Q. Ma, I. Schmidt, J. Soffer, and J.- J. Yang, Phys. Lett. B 547, 245 (2002). 



16 




FIG. 11: (Color online) The ratio /p^ (right) and p^/p'^(left) at = 5 GeV^ with q = 
0.25 GeV/fm^. The solid, dashed and dotted curves are the predictions of the cluster model, 
the rescaling model and the pion excess model respectively. Thin, normal and thick curves corre- 
spond to the three options of the unfavored fragmentation (1), (2), (3) respectively. The target 
nuclei assumed here is Fe. 

[24] B.-Q. Ma, I. Schmidt, and J.-J. Yang, Phys. Lett. B 477, 107 (2000). 

[25] A. Bialas and T. Chmaj, Phys. Lett. B 133, 241 (1983); A. Bialas and M. Gyulassy, Nucl. 

Phys. B291, 793 (1987). 
[26] E. Wang and X.N. Wang, Phys. Rev. Lett. 89, 162301 (2002). 
[27] X. Guo and X.N. Wang, Phys. Rev. Lett. 85, 3591 (2000); 

X.N. Wang and X. Guo, Nucl. Phys. A 696, 788 (2001); 

J. A. Osborne, E. Wang, and X.N. Wang, Phys. Rev. D 67, 094022 (2003). 
[28] B. Kopehovich, J. Nemchik and E. Predazzi, | n"ucl-th/96070"36 
[29] A. Accardi and H.J. Pirner, Nucl. Phys. A 711, 264 (2002). 
[30] X.N. Wang, Z. Huang, and 1. Sarcevic, Phys. Rev. Lett. 77, 231 (1996). 
[31] F. Arleo, Eur. Phys. J. C 30, 213 (2003). 
[32] F. Arleo, JHEP 0211, 044 (2002). 



17 




FIG. 12: (Color online) The ratio {p^ / p^) / {p^ / p^) at = 5 GeV^ with q = 0.72 GeV/fm^ (left) 
and q = 0.25 GeV/fm^ (right). The solid, dashed and dotted curves are the predictions of the 
cluster model, the rescaling model and the pion excess model respectively. Thin, normal and thick 
curves correspond to the three options of the unfavored fragmentation (1), (2), (3) respectively. 
The target nuclei assumed here is Fe. 
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FIG. 13: (Color online) The ratio A^/A^ at = 5 QeV^ with q = 0.25 GeV/W. The solid, 
dashed and dotted curves are the predictions of the cluster model, the rescaling model and the pion 
excess model respectively. The left figure is the result in the region 0.3 < z < 0.5 and the right is 
the result in the region 0.7 < z < 0.9, with the thin, normal and thick curves corresponding to the 
three options of the unfavored fragmentation (1), (2), (3) respectively. The target nuclei assumed 
here is Fe. 
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FIG. 14: (Color online)The production ratio of A^/A^ at = 5 GeV^ with q = 0.25 GeV/W. 
The sohd, dashed and dotted curves are the predictions of the cluster model, the rescaling model 
and the pion excess model respectively. The left figure is the result in the region 0.3 < z < 0.5, 
the middle is the result in the region 0.4 < z < 0.9 and the right is the result in the region 
0.7 < z < 0.9. The thin, normal and thick curves correspond to the three options of the unfavored 
fragmentation (1), (2), (3) respectively. The target nuclei assumed here is Fe. 
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FIG. 15: (Color onlme)The production ratio of (A^/A^)/(A^/A^) at = 5 QeV^ with q = 
0.25 GeV/fm^, calculated in 0.7 < z < 1.0 and 0.4 < z < 1.0. The solid, dashed and dotted curves 
are the predictions of the cluster model, the rescaling model and the pion excess model respectively. 
The thin, normal and thick curves correspond to the three options of the fragmentation (1), (2), 
(3) respectively. The target nuclei assumed here is Fe. 
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